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ABSTRACT 

We present number counts, luminosity functions (LFs) and luminosity densities 
of galaxies obtained using the Sloan Digital Sky Survey Sixth Data Release in all 
SDSS photometric bands. Thanks to the SDSS DR6, galaxy statistics have increased 
by a factor of ~ 9 in the M-band and by a factor of ~ 4 — 5 in the rest of the SDSS 
bands with respect to the previous work of lBlanton et alj (|2003bD . In addition, we have 
achieved a high redshift completeness in our galaxy samples. Firstly, by making use 
of the survey masks, provided by the NYU-VAGC DR6, we have been able to define 
an area on the sky of high angular redshift completeness. Secondly, we guarantee that 
brightness-dependent redshift incompleteness is small within the magnitude ranges 
that define our galaxy samples. With these advances, we have estimated very accurate 
SDSS DR6 LFs in both the bright and the faint end. In the 01 r-band, our SDSS DR6 
luminosity function is well fitted by a Schechter LF with parameters $* = 0.90 ±0.07, 
M*—5logioh = — 20.73±0.04 and a = — 1.23±0.02. As compared with previous results, 
we find some notable differences. In the bright end of the 01 u-band luminosity function 
we find a remarkable excess, of ~ 1.7 dex at Mo.i„ ~ —20.5, with respect to the best- 
fit Schechter LF. This excess weakens in the 01 <7-band, fading away towards the very 
red 2-band. A preliminary analysis on the nature of this bright-end bump reveals 
that it is mostly comprised of active galaxies and QSOs. It seems, therefore, that 
an important fraction of this exceeding luminosity may come from nuclear activity. 
In the faint end of the SDSS DR6 luminosity functions, where we can reach 1 — 1.5 
magnitudes deeper than the previous SDSS LF estimation, we obtain a steeper slope, 
that increases from the 01 u-band, with a = —1.01 ± 0.03, to the very red 01 z-band, 
with a = — 1.26±0.03. These state-of-the-art results may be used to constrain a variety 
of aspects of star formation histories and/or feed-back processes in galaxy formation 
models. 

Key words: catalogues - surveys - galaxies: luminosity function, mass function - 
large-scale structure of Universe. 



1 INTRODUCTION 



From the pioneering work of iHumasonl (|l956l ) and ISandagel 
1 19781 ). who measured redshifts of bright galaxies from th e 
Shapley-Ames photometric catalog Tshaplev fc Amesl l|l932l ). 
much effort has been invested in mapping the luminous 
and matter contents of the Universe. The Center for As- 
trophysics survey (CfA, iHuchra et alj Il98&t ) is considered 
as the first proper redshift survey, specifically designed for 
cosmology studies. More important, from the first slice of 
about 1,000 galaxies the CfA Redshift Survey provided 
the community with observational evidence of an old the- 
oretical, and at times controversial idea: the e xistence of a 
large scale structure of galaxies in the Universe (iDavis et al.l 
1 19851 ). This first vision of cosmic complexity encouraged 
the development of new imaging and spectrometric technol- 
ogy and, consequently, gave rise to a number of other red- 



shift surveys that followed different approaches and strate- 
gies. To na me but a few, the So uthern Sky Redshift Sur- 
vey ( SSRS, da Costa et al. 1988); the Perseus-Pisces cat- 
alog (|Giovanelli fc Havnesl Il99ll ) or the catalogs based on 
data from the Infrared Astronomical Satellite (IRAS). In 
the last decade, the emergence of multi-fiber spectrographs 
set the scene for larger and deeper redshift surveys. Ex- 
amples of these are the Las Campanas Redshift Survey 
(LCRS, IShectman et al.lll99rj ). consisting of 26,418 galax- 
ies with an average redshift of z ~ 0.1; or the 2 degree 
Field Redshift Survey (2dF RS, IColless et al.ll200ll ). with 
about 222, 000 galaxies and covering a sky area o f 1500 deg 2 . 
Final ly, the Sloan Digital Sky Survey (SDSS, lYork et all 
2000) is the largest photometric and spectroscopic survey 
ever compiled, and represents the most accurate map of the 
nearby universe at z < 0.3. The SDSS Sixth Data Release 
ijAdelman-McCarthv et al]|2008l ). that we use in this paper, 
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contains spectroscopic information for more than 1,000,000 
galaxies and quasars which spread over 7425 deg 2 on the 
sky. Only in recent yea rs, with surveys like the DEEP2 
Galaxy Redshift Survey (iDavis et al.ll2003l) or th e VIMOS- 
VLT Deep Survey fVVDS. iLe Fevre et al.l l2003h . have we 
reached the stage where we can study the galaxy popula- 
tion in the distant (z ~ 1) Universe. Other high-z surveys 
are currently being completed. 

The advances in the survey field also made it necessary 
to develop data reduction pipelines and analysis tools to 
process and understand increasingly larger data sets. These 
days, cosmologists use a number of statistics to characterize, 
for a particular survey, the distribution of galaxies in three- 
dimensional space. Number counts, selection functions, lu- 
minosity functions or correlation functions are just a few 
examples. In this work we focus on the number counts and 
the luminosity functions of galaxies, that we draw from the 
SDSS DR6. Number counts, which describe the distribution 
of fluxes of galaxies, have been calculated for a number of 
surveys. The general consensus is that, in the close- by uni- 
verse, galaxy number counts look l ike what we expect from 
an euclidean, not-evolved Universe. lYasuda et al.l l|200ll ) ob- 
tained number cou nts for the SDSS Co mmissioning Data 
in all ugriz bands. iNorberg et ail (2002) pro vided number 
count s for the 2dF survey in the bj band. iFeulner et al.l 
(2007) not only estimated galaxy number counts for a set of 
catalogs b ased on the Munich Near-Infrared Cluster Survey 
CMUNICS. lDrorv et alllioOlh but also presented a complete 
revision of this subject in the literature (see their Figure 8). 
In contrast to galaxy number counts, the luminosity func- 
tion (LF), which is the number density of galaxies per unit 
absolute magnitude, ha s been historically a rather controver- 
sial is sue. For examp le, Marzke et al.l (|l994T ) - using the CfA 
RS -. INorberg et all {{200$ ) - 2dF RS - and iBlanton et ail 
(|2003bl ) - SDSS DR2-. all obtained very different results. 

Both the luminosity function and the luminosity density 
of galaxies are observational signs of the process of galaxy 
formation and evolution. A precise determination of these 
statistics is needed to constrain current theories. Conse- 
quently, new discoveries in observational cosmology could 
make a strong impact in our understanding of the physi- 
cal processes that drive the birth and life of galaxies in the 
Universe. Nowadays, state-of-the-art models of galaxy for- 
mation invoke a number of galactic "mechanisms", which 
are connected through the so-call feed-back processes. Dis- 
entangling these relations is an ambitious but crucial task in 
modern Cosmology. In this sens e, the semi-analytic models 
of galaxy formation (SAMs, e.g. ICroton et al]|2006l ). which 
are emb edded in N-body sim ulations like the Millennium 
Run (see lSpringel et al.ll2005l ). are a very useful tool for cos- 
mologists. These SAMs are a good ground for testing new 
theoretical ideas and understanding their observational im- 
plications. 

The main purpose of this work is to take advantage 
of the large increase in the galaxy statistics thanks to the 
SDSS DR6 to obtain the number counts, LFs and luminos- 
ity densities of galaxies in the nearby universe. We intend 
to shed light both in the faint end of the LF, where most 
discrepancies come from, and in the bright-end, where statis- 
tics have always been poor and errors, consequently large. 
In section [2] we briefly describe the SDSS DR6, discuss our 
sample selection and comment on redshift completeness. In 



section [3] we present our results on the number counts, the 
LFs and the luminosity densities of galaxies in each one of 
the SDSS photometric bands. Finally, in Section [3] we dis- 
cuss our results and in Section [5] we present a summary of 
our work. Throughout this paper, unless otherwise stated, 
we assume a standard ACDM concordance cosmology, with 
Sl m = 0.3, Qa ~ 0.7, w = —1, and h = 1. In addition, we 
use AB magnitudes. 



2 THE SLOAN DIGITAL SKY SURVEY DATA 
RELEASE 6. DATA SAMPLES SELECTION 
AND REDSHIFT COMPLETENESS 

In this work, we use the SD SS Sixth Data Release 
l|Adelman-M cCarth v et al"1l2008h . This data set completes 
the North Galactic Cap, containing photometric informa- 
tion of ~ 290 million objects over 9583 deg 2 . Around 1.27 
million objects were selected for spectroscopy, covering an 
area of 7425 deg 2 on the sky. Important for this work, spec- 
troscopy is availab le for ~ 1, 000, 000 galaxies down to mag - 
nitude r ~ 17.77 (|York et al.ll2000l ; IStoughton et al.ll2002r i. 
Detailed information about the SDSS DR6 can be found in 
lAdelman- McCarthy et all (120081 ). 

The SDSS DR6 is the largest survey of the nearby 
universe publicly available. The SDSS collaboration have 
successively extended their catalogs since the times of the 
SDSS Early Data Release and, consequently, improved enor- 
mously our capability of ma pping the univ e rse up to redshift 
2 ~ 0.25. A few years ago, IBlanton et alj (|2003bl ) used the 
SDSS DR2 to estimate the luminosity function of galaxies. 
In this work, thanks to the SDSS DR6, the size of our sam- 
ples has risen by a factor ~ 9 in the very blue u-band and 
by a factor ~ 4 — 5 in the other bands ( g, r, i and z). This 
huge enhancement in the statistics will be especially critical 
in the bright end of the LF, where statistics have always 
been poor in the past. In addition, we also expect to reach 
deeper magnitudes in the faint end with respect to previous 
works. Galaxy number counts can also be estimated with 
significantly more accuracy. It is therefore well justified to 
update the current knowledge on the number counts, lumi- 
nosity functions and luminosity densities of galaxies in the 
close-by Universe. 



2.1 Data samples selection 

In this section, we describe the selection of our galaxy sam- 
ples. We have drawn our samples from the NYU Value 
Added Galaxy Cata log DR6 Large Scale Structure sample 
l|Blanton et al.ll2005l ). The NYU-VAGC is a compilation of 
galaxy catalogs cross-matched to the SDSS that includes a 
number of useful quantities derived from the photometric 
and the spectroscopic catalogs (such as K-corrections or ab- 
solute magnitudes). It also incorporates a precise and user- 
friendly description of the geometry of the survey, which has 
never been used in previous SDSS LF works. 

Before we proceed with the sample selection, we will 
describe how we built our Parent Sample from the NYU- 
VAGC LSS sample. The Parent Sample is the galaxy cat- 
alog from which we extract all samples that we use in this 
work. Firstly, we take all objects that satisfy the criteria of 



The SDSS DR6 Luminosity Functions of Galaxies 3 




the SDSS Main Galaxy Sample l|Strauss et alj|2002l ). Sec- 
ondly, we make use of the survey masks, provided by the 
NYU-VAGC, to define our high- completeness survey win- 
dow (hereafter, HCSW). In this paper we deal with number 
counts and luminosity functions in the entire SDSS DR6 
spectroscopic catalog. These tasks require knowing the ef- 
fective area covered on the sky by the sample . In practice, 
redshift completeness, which is the fraction of galaxies with 
a successful redshift measurement, is far from being uniform 
across the sky. In Figure [1] we show the angular sky redshift 
completeness of the SDSS DR6 spectroscopic catalog for a 
major part of the survey (~ 99%), just excluding objects 
outside the range 110° < RA < 270°, for the sake of clarity. 
The plot has been pixelized so that each pixel - or poly- 
gon - is an area of constant completeness. Polygons are also 
color-coded in white to black tones, the latter meaning 100% 
completeness. In order to avoid overestimating the area cov- 
ered on the sky by our samples, the HCSW is defined so that 
redshift completeness is at least 80% in every polygon. Only 
galaxies lying in these polygons are included in our Parent 
Sample. 

With these first restrictions our Parent Sample is com- 
prised of 947, 053 galaxies that spread over 6428.32 deg 2 
on the sky. From this catalog, each sample is drawn by 
applying the following cuts to the redshift and apparent 
magnitude: 

• n% m i n 

(j) < m CD < mmaxiS) 

• Z m in($) < z(j) < ZmaxQ) 

where j=u, g, r, i, z. In Table [1] we show lower and upper 



Band 


Number 


'fTlmin 




z rnin 


Zmax 


u 


192,068 


16.45 


19.00 


0.02 


0.19 


g 


241,719 


14.55 


17.91 


0.02 


0.16 


r 


516,891 


13.93 


17.77 


0.02 


0.22 


i 


429,173 


13.55 


17.24 


0.02 


0.22 


z 


414,828 


13.40 


16.97 


0.02 


0.23 



Table 1. Number of galaxies and limits in apparent magnitude 
and redshift of each SDSS photometric band sample. Motivation 
for each cut is discussed in section [2. II 

limits of these quantities along with the number of galaxies 
for each sample. 

At this point, it is convenient to clarify the motiva- 
tion for each cut. The apparent magnitude limits of Ta- 
ble [T] are set to ensure that the effect of redshift incom- 
pleteness is small in our SDSS galaxy samples. In the bright 
end, brightness-dependent redshift incompleteness starts to 
be important at r < 15. In the faint end, redshift incom- 
pleteness in all SDSS bands is domin ated by the intrinsi c 
faint limit of the Main Galaxy Sample (|Strauss et alj|2002h , 
i.e. r = 17.77. In order to choose these apparent magni- 
tude limits, we have made use of the SDSS galaxy num- 
ber counts, that will be properly discussed in Section 13.11 
For each galaxy sample, we have taken the magnitude range 
where the number of galaxies rises at constant rate in each 
SDSS band, with a 0.1-dex deviation allowance (see Figure 
4, where galaxy number counts have been scaled by an eu- 
clidean, not-evolved model). Within these limits, we have 
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Figure 2. Redshift distributions in all SDSS photometric bands. 
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Figure 3. Redshift completeness vs. apparent magnitude in the 
r-band. At r < 14 completeness decreases sharply. Errors have 
been estimated by propagating the poissonian uncertainties to 
the redshift completeness. 



data i s estimated in a 6% incompleteness bv lStrauss et al.l 
(|2002h . 

In Figure [3] we show how redshift completeness varies 
with apparent magnitude in the r-band. At r ~ 14, redshift 
completeness falls from ~ 85% to ~ 50% at r ~ 12. In the 
faint end, however, it reaches a plate au at 90% - down t o 
r ~ 17.8 -, which is in agreement with lStrauss et all l)2002h . 
In the rest of the SDSS bands, redshift completeness not only 
decreases in the bright end, but also in the faint end. This 
is due to the intrinsic r-band faint limit of the Main Galaxy 
Sample at r = 17.77. Interestingly, the apparent magnitude, 
m-min(j), at which this decrease occurs varies between bands 
as a result of the dispersion in the colors of galaxies. We have 
checked that, by imposing m max (j) > m(j) > m mi „(j), 
we ensure that each galaxy sample is approximately 85% 
complete in any magnitude bin. 

Another aspect of redshift completeness that should be 
carefully taken into account is its angular variation. As we 
discussed above, a precise estimation of the area covered in 
the sky by our samples is needed in order to calculate num- 
ber counts and luminosity functions of galaxies. We ensure 
here that we do not overestimate this area by defining a re- 
gion, the HCSW, in which redshift completeness is at least 
80%. The HCSW covers 6428.32 deg 2 , about 90% of the en- 
tire sky coverage of the SDSS DR6 spectroscopic catalog. 
As we can see in Figure [T] the SDSS DR6 sky coverage is 
patched with areas where no redshift was measured at all. 
These regions cover a total of 623.5 deg 2 on the sky. 

Finally, we cannot discard the possibility that the SDSS 
spectroscopic catalog is incomplete f or very low surface 
brightness objects, i.e. fi r ,50 > 24. fsee lStrauss et aT] I2002T ) 
The presence of this selection effect in the data could in prin- 
ciple affect our results. However, we have evidence that the 
surface brightness of most galaxies in the SDSS, and conse- 
quently in our galaxy samples, is far greater than (i r ,5o — 24 
|Blanton et alj|2003b| ). 



estimated that redshift completeness is ~ 85% in all SDSS 
bands. In Section ^. 21 we will discuss on redshift incomplete- 
ness issues in the SDSS in more detail. 

In Figure [2] we show the redshift distribution for the 
Parent Sample in each SDSS photometric band. This figure 
illustrates the motivation for the redshift limits given in Ta- 
ble [1] The lower redshift limit is set to z = 0.02 to avoid 
the redshift incompleteness that affects the very bright and 
nearby galaxies. The upper redshift limit corresponds to the 
redshift at which 95% of objects are selected in each sample 
and is set for consistency. 

2.2 Redshift Completeness 

It is well known that the brightest galaxies in the SDSS 
catalogs are affected by a severe redshift incompleteness. 
This effect is especially important for nearby galaxies. Ap- 
parently large and complex objects represent a major task 
for the photometric SDSS pipelines l|Strauss et alj|2002t ). In 
addition, and affecting the entire magnitude range, is the 
so-called fiber collisions problem. This source of incomplete- 
ness is due to the fact that fibers cannot be placed closer 
than 55". The consequence of this limitation in the SDSS 



3 RESULTS 

3.1 Number counts 

In Figure [4] we plot with different symbols the logarithm 
of the number of galaxies per unit area and apparent mag- 
nitude (actually, half magnitude), scaled by an euclidean 
model for all SDSS bands. We choose, arbitrarily, the fol- 
lowing model for euclidean counts: 

N euclidean = io°- 6 <»- 18 > (1) 

where x is replaced by the apparent magnitude in each SDSS 
band. 

In the r-band, galaxy number counts increase by a factor 
of about 10 from r ~ 12 to r ~ 13.5. This magnitude range 
is strongly affected by redshift incompleteness, as discussed 
above (see also Figure|3]). From r ~ 14 to r ~ 18, counts rise 
at approximately the same rate as that of the model (Equa- 
tion [Q. At r ~ 18, where the SDSS spectroscopic faint-end 
limit is set (r = 17.77), galaxy number counts fall sharply. 

In the rest of the bands, the behavior is very similar in 
the bright end and it is also due to redshift incompleteness. 
In the z-band, counts start to follow the euclidean model at 
z ~ 13.5. In the i-band, this happens at i ~ 13.5; in the 
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Band 


Apparent magnitude limits 


Rcdshift limits 


Slope 


u 


16.45 < u < 19.00 


0.02 < z < 0.19 


0.627 


g 


14.55 < g < 17.91 


0.02 < z < 0.16 


0.608 


r 


13.93 < r < 17.77 


0.02 < z < 0.22 


0.610 


i 


13.55 < i < 17.24 


0.02 < z < 0.22 


0.611 


z 


13.40 < z < 16.97 


0.02 < z < 0.23 


0.626 


Table 


2. Slope in units of mag -1 


of galaxy number 


counts 



within the apparent magnitude and redshift limits of each sam- 
ple. Within these ranges, number counts are consistent with an 
Euclidean, not-evolved Universe. 
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Figure 4. Galaxy number counts in all SDSS bands scaled by 
an euclidean model in bins of half a magnitude. Poissonian errors 
are of similar size as symbols, so they are not shown. 



g-band, at g ~ 14.5 and, in the u-band, at u ~ 16.5. In 
the faint end, galaxy number counts fall at u ~ 19, g ~ 18, 
i ~ 17 and z ~ 17. This decrease is obviously less pro- 
nounced than in the r-band. This reflects the fact that the 
color-magnitude relation is not univoque. In addition, the 
dispersion in this relation is con siderably larger for u — r 
colors than it is for r — z colors (jBlanton et al.ll2003al ). As 
seen in Figure [4] the slope of the faint-end decrease in the 
galaxy number counts is considerably steeper in the u and 
g-bands than it is in the redder bands. 

Number counts are consistent with an Euclidean, not- 
evolved Universe in all SDSS bands within the magnitude 
and redshift ranges given in Table [1] Within these ranges, 
the slopes in units of mag -1 of the SDSS DR6 galaxy num- 
ber counts are listed in Table [2] Although discrepancies ex- 
ist with respect to the euclidean model , these are probably 
smaller than the expected uncertainty in the determination 
of the apparent magnitude and redshift of each galaxy in 
the SDSS. 

Our results for the nearby univers e are in agreement 
with a number of previous works (see Feulner et alj 120071 
for a review). However, only lYasuda et al.l <|200lf > obtained 
galaxy number counts for SDSS data. The authors used 
imaging data taken only during the commissioning phase of 
the SDSS. They found an euclidean- like behavior up to mag- 
nitude m min ~ 12 - except for the u-band, where m min ~ 14 
- (see their Figure 8) . The lack of galaxies in the very bright 



end of Figure [4] with resp ect to photometric number counts 
from lYasuda et all ([2001) is partially due to the strong red- 
shift incompleteness that affects the SDSS spectroscopic cat- 
alog in these magnitude ra nges (see Sect i on 12.211 . However, 
our results and those from lYasuda et ahl |200lf) are not di- 
rectly comparable. Firstly, they used a very limited sample, 
in terms of sky coverage (~ 230 deg 2 ). Secondly, comparing 
spectroscopic and photometric results is always tricky. 



3.2 Luminosity Functions 

In order to estimate the luminosity function of galaxies in 
each SDSS photometric band, we take absolute magnitudes 
and K-correct ions from the NYU V AGC DR6 LSS cata- 
log. Following iBlanton et ail (|2003bh . absolute magnitudes 
are calculated with the SDSS photometric bands shifted to 
z = 0.1. With this convention, the absolute magnitude of a 
galaxy in a given band shifted to z = 0.1 
constructed from its apparent magnitude at z 
its redshift z as follows: 



Mo.ij, would be 
0, rrij, and 



Mo, 



5log 10 h - DM(z) - Ka,ij(z) 



(2) 



where DM(z) is the distance modulus (which de- 
pends also on the cosmological parameters) and Ko,xAz), 
the K-correction for the galaxy in the shifted band 01 j. 
IBlanton et aD l|2003bh included another correction in Ex- 
pression [2] to account for the evolution of the luminosity of 
a galaxy from redshift z to z = 0, the so-called evolution cor- 
rection. However, we expect this correction to be very small 
within the redshift ranges that we have considered here (see 
also Section [4| . In Figure [5] we show the distribution of K- 
corrected absolute magnitudes in each galaxy sample. The 
shape of the absolute magnitude distribution is very simi- 
lar in all SDSS bands: a gaussian-like distribution slightly 
skewed to fainter magnitudes. However, mean values move 
towards bright-end bins from 0,1 it-band - Mo.i„ ~ —18 - to 
01 z-band - Mo.i z ~ —21.5, which is consistent with the fact 
that red objects are, on average, brighter than blue objects. 
In Figure [6] we also show the bimodal 01 u — 01 r color dis- 
tribution of galaxies in our 01 r-band sample. With a dashed 
line we represent the demarc ation commonly used to sepa- 
rate red and blue objects (see IStrateva et a"fll200lf ). 

We use the Stepwise Ma ximum Likelihood Method 
(SWML: lEfstathiou et al.lll988h to estimate the LF of galax- 
ies, that is commonly expressed as &(L). This technique does 
not rely on any assumption about the shape of $(!/) and 
is considered as the most reliable e stimator of the LF. For 
mor e details about this m ethod see lEfstathiou et alj (|l98&t ) 
and lNorberg et al.l (|2002l ). The SWML requires an indepen- 
dent estimation of the number density of galaxies, i.e. the 
normalizatio n constant, n. We use th e following prescription 
proposed by iDavis fc Huchral l| 19821 ) based on the selection 
function of each galaxy, 4>{ z i)y an d the maximum volume 
encompassed by the sample, V m ax- 



E; 



(3) 



In Figure [7] we show the SWML estimate of the SDSS 
DR6 LF in the 01 r-band. In addition, we over-plot with a 
dashed line the best-fit Schechter function ilSchechtedli976h . 
which has the following shape: 
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Figure 5. Absolute magnitude distribution in each sample. Al- 
though the shape of the distributions is very similar, mean values 
move towards bright-end magnitude bins from 01 iz-band to 01 2- 
band, as red objects are, on average, brighter than blue objects. 
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Figure 6. The ( 01 «— 01 r) color-magnitude diagram in the 01 r- 
band. The dashed line represents the demarcation commonly used 
to separate red and blue galaxies. 



$(M) = o.4io g (io)$ao-°' 4(M - M * )(Q+1) 

expf-lO- ' 4 ^-"*'] (4) 

where a, M* and $* are the three parameters to fit. 
Values of these parameters for the best-fit Schechter func- 
tion are given in Table [3] For comparison , we also show in 
Figure \7\ the LF of iBlanton et al.l (|2003bl ^ with a solid line. 
This comparison will be addressed in the Discussion sec- 
tion. To calculate errors in the SWML estimates of the LF 
we perform a bootstrapping analysis using f,000 random 
sub-samples of 1/3 of the number of objects in each sample. 
In Figure [7] shaded regions represent the la uncertainty ob- 
tained from this method. 

Due to the big number statistics that we have, with 
about 500, 000 galaxies in the 01 r-band, errors are only sig- 
nificant in the very bright end of the LF. In the faint end, 
we can go down to Mo.i r ~ —16, which means that we can 
build the LF with unprecedented precision within a very 



Figure 7. The 01 r-band SDSS DR6 Luminosity Function. The 
SWML LF estimate is shown in diamonds. The dashed line rep- 
resents the best-fit Schechter function and the solid line, the 01 r- 
band LF from Blanton et al. (2003b). Best-fit values of Schechter 
parameters a, M* and 3> t are also shown in the figure. Shaded 
regions represent the lcr uncertainty calculated using a bootstrap- 
ping technique. 

large range of magnitudes. As we will see below, the above 
statements hold for all SDSS bands. Our 01 r-band LF is 
reasonably well fitted by a Schechter LF with a faint-end 
slope a — —1.23. It is only in the very bright end where 
this best-fit Schechter LF starts to underestimate our LF. 
At Mo.i r < —23.5, statistics are poor and errors become 
increasingly large. 

In Figure [8] we present, in the same way as in Figure [TJ 
SWML estimates of the LF in bands V l g, aA i and °- l z, 
as well as their corresponding best-fit Schechter LF. Values 
of best- fit Schechter parameters are also given in Table [3] As 
in the r-band, errors are only significant in the very bright 
end of °' 1 u, g, i and 0,1 z - band luminosity functions. In 
addition, we can go down to very faint magnitudes without 
losing precision. 

In the very blue 01 u-band, the shape of our SDSS DR6 
LF is consistent with a Schechter LF with an almost flat 
faint-end slope (corresponding to a = —1.01). However, in 
the bright end, we find a remarkable excess with respect 
to the best-fit Schechter LF. This excess, of ~ 1.7 dex at 
Mo.i" ~ —20.5, is very significant within the magnitude 
range -20.5 < Mo.i„ < -22. In the 01 g-band, the bright- 
end bump weakens considerably, but it is probably still sig- 
nificant, even though errors are large according to our boot- 
strapping analysis. In this band, our SDSS LF is very well 
fitted by a Schechter LF with a positive faint-end slope, cor- 
responding to q = —1.10. Only in the very bright-end, where 
the excess is still noticeable, do we find some discrepancy. 
Below, we provide a preliminary analysis and discussion on 
the nature of the 01 u-band LF bright-end bump, that may 
have important implications in terms of galaxy formation 
and evolution. 

In the redder bands we find a positive faint-end slope, 
corresponding to a = —1.16 in the 01 i-band and a — —1.26 
in the 01 ,z-band. The bright-end bump has diminished sub- 
stantially in the 01 i-band and has disappeared in the very 
red 01 z-band. It is interesting to note that, from the u- 
band to the 01 z-band, the shape of the SWML estimate of 
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Figure 8. The SDSS DR6 Luminosity Functions in bands 01 «, 01 <;, 01 i and 01 2. SWML LF estimates are shown in diamonds and 
best-fit Schechter functions are represented by dashed lines. In addition, we over-plot in each panel the LF from Blanton et al. (2003b). 
Best-fit values of Schechter parameters a, M* and <£>» are also shown. Shaded regions represent the la uncertainty calculated using a 
bootstrapping technique. 



the SDSS LF changes following a clear pattern. The faint- 
end slope increases towards the redder bands (see Table [3}, 
being almost flat in the 01 u-band and remarkably steep in 
the 01 2-band. Conversely, the excess of bright galaxies with 
respect to the best-fit Schechter function diminishes towards 
redder bands. In principle, one might think that the strong 
bright-end excess in the 01 «-band fades progressively to- 
wards redder bands and only disappear completely in the 
z-band. In this sense, the ai r-band SDSS LF seems to 
slightly deviate from this trend. In the faint-end, we find a 
slope that is a bit larger than we could expect (a — —1.23), 
but this could be just a consequence of the fact that the 
entire SDSS spectroscopic sample was selected in this band. 

The bright-end bump that shows up clearly in the 01 it- 
band LF (and partially in the 01 g-band LF ) is an interesting 
discovery that may have implications for our understanding 
of galaxy formation and evolution. In order to investigate 



the nature of the objects that populate it, we have selected 
all galaxies brighter than —20.5 in the 01 w-band sample. We 
will hereafter refer to this sample, which is comprised of 422 
objects, as the bright-end bump sample (BEBS). 

We have visually inspected the spectra of all galaxies 
in the BEBS and classified them into three types. About 
60% of objects have a typical QSO or Seyfert I spec- 
trum, 27% of sources are classified as LINERs or Star- 
burst (SB) and only ~ 13% of objects are galaxies with 
a bulge-like spectrum, and showing little or no signs of 
star formation. In Figure [9] we plot in the (°" u — °" r) 
vs. Mo.i u colour-magnitude diagram (CMD) of the °' 1 u- 
band sample the three types of BEBS galaxies discussed 
above. In general, these galaxies form a relatively tight se- 
quence in the bright-end of the CMD, showing a consid- 
erable color dispersion. Note that, in contrast to Figure |6j 
density contours are now log-spaced and hence, these ob- 
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Table 3. Values of Schechter parameters a, M* and <I>* of the best-fit Schechter function in all SDSS bands for this work and for Blanton 
et al. (2003b). 



jects occupy an extremely underpopulated region in the 
CMD. In the left-hand plot, QSOs/Seyferts I show the 
smallest color dispersion and are, on average, the bluest: 
< ( 01 u - 01 r) >= 0.68. Both LINERs/SBs (middle plot) 
and bulge-like galaxies (right-hand plot) show much larger 
color dispersion and are, on average, considerably redder, 
with < { 0A u - 01 r) >= 1.07 and < ( 01 u - 01 r) >= 0.91, 
respectively. In the 01 u-band, bulge-like galaxies are, on av- 
erage, the brightest among the BEBS galaxies: < Mo.i u >= 
(-20.97 QSO/S/ , -20.77 i/JV/ss , -21.03 bul3e ). 

We have checked that BEBS galaxies are typically at 
high redshift relative to the average redshift of the 0,1 it-band 
sample from which they are drawn. The mean redshift in the 
BEBS is < z >~ 0.155 while in the entire 01 u-band sample 
is < z >~ 0.080. 

At this point, it is necessary to remark that this is a 
preliminary analysis. We have performed a rough classifica- 
tion of the BEBS galaxies, and this is still subject to some 
uncertainty. However, our aim is to provide a first approach 
to the nature of this population of objects. In Section [4] we 
speculate on possible implications of these results. 

Finally, in Figure 1101 we present, in the same way as 
in Figure \7\ and Figure O SWML estimates of the SDSS 
0,1 r-band LF for blue and red galaxies separately, as well 
as their corresponding best-fit Schechter LF. Values of best- 
fit Schechter parameters are also shown in the figure. The 
demarcation that we use to separate blue and red objects 
is represented by a dashed line in Figure [6] The SDSS r- 
band LF of blue galaxies is well fitted by a Schechter LF 
with a = -1.36. In contrast, the SDSS Vband LF of red 
galaxies has a negative faint-end slope, corresponding with 
a = —0.81. These results are consistent with a number of 
previous works. In the bright end, at Mo.i r < —21.5, the 
blue LF falls remarkably below the red LF. 



3.3 Luminosity Densities 

To estimate the luminosity density, p, in all SDSS bands, we 
integrate throughout the entire absolute magnitude ranges 
defined by the limits shown in Table [4] In this table we list 
p in all bands obtained using the SWML estimate of the 
luminosity functions shown in Figure [7] and Figure [5] Lu- 
minosity densities are expressed in AB mags 
We also give in Table|4]the effective wavelength correspond- 
ing to each SDSS photometric band, \ e ff- Errors shown in 
this table have been calculated using a similar bootstrap- 
ping technique as that discussed in Section 13.21 We have 



used 1000 random sub-samples of 1/3 the number of objects 
in each sample. 

As expected, the luminosity density in the nearby uni- 
verse increases with \ e ff , in absolute values. This means 
that the Universe, up to z ~ 0.2, is considerably more lumi- 
nous in the red side of the spectrum than it is in the blue 
side. Our luminosity densities obtain ed using the SWML LF 
estimate are in good agreement with iBlanton et all l|2003bl ) 
(see Table |3| . Both estimates differ in less than 2% in any 
band. 

However, it is necessary to note that the value of p is 
not too sensitive to small variations in the shape of the lu- 
minosity function. In Table 3J we also provide luminosity 
densities in each band obtained using the Schechter fits to 
the SWML LF estimates of the luminosity functions dis- 
cussed in Section 13.21 The values of p obtained using the 
SWML estimates and those obtained from the Schechter fit 
differ in lees than 1% in each band. 



4 DISCUSSION 

The main results presented in this work are the SDSS DR6 
Luminosity Fu nctions of galaxies in the nearby universe. A 
few years ago, IBlanton et al.l l|2003bl ) used an early version 
of the SDSS (DR2) to calculate the SDSS galaxy LFs. Now, 
with the SDSS DR6 available, galaxy statistics have im- 
proved by a huge factor of ~ 9 in the very blue 01 w-band and 
by a factor of ~ 4 — 5 in the rest of the SDSS photometric 
bands. Moreover, we have ensured a high redshift complete- 
ness in our galaxy samples. Firstly, we have defined an area 
on the sky of high angular redshift completeness by making 
use of the survey masks, provided by the NYU-VAGC. Sec- 
ondly, we guarantee that the effect of brightness-dependent 
redshift incompleteness is negligible within the magnitude 
ranges that define our galaxy samples. These advances make 
our SDSS DR6 LF s substa ntially more precise than those 
from IBlanton et all (|2003bT ) at both the bright and the faint 
end. This said, the LFs of IBlanton et all l|2003bl ) are com- 
patible with our results. However, notable differences, which 
are surely physically significant, exist. In the bright end of 
the blue bands LFs (especially in the 0,1 it-b and) we find a 
remark able excess, which was very noisy in IBlanton et al.l 
l|2003bl ) due to their lack of statistics. In the faint end, we 
obtain steeper slopes in all SDSS bands, especially in the 
01 M-band, where the DR6 statistics allow u s to go about 
1 — 1.5 magnitudes deeper as compared to IBlanton et al.l 
(2003b). If the huge improvement in the statistics and/or 
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Table 4. Luminosity densities in AB mags (Mpc^h^ 1 )^ 1 for all SDSS photometric bands calculated using the SWML estimate of the 
luminosity functions of Figure [7] and Figure [8] (pswml) and the Schechter LF best fits of Table \3\ (ps'chechter)- Absolute magnitude 
ranges of the integration and the effective wavelength corresponding to each band, \ e j j , are also provided. Errors in Pswml have been 
calculated using a bootstrapping technique. In addition, we give for comparison the luminosity densities from Blanton et al. (2003b). 
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Figure 9. The ( 0A u — 01 r) vs. Mo.i u color-magnitude dia- 
gram for the three types of bright-end bump galaxies considered: 
Seyferts I/QSO's (left-hand plot), LINER's/SB's (middle plot) 
and Bulge-like galaxies (right-hand plot). The underlying CMD 
of the entire 01 u-band sample from which BEBS galaxies are 
selected is shown with log-spaced contours. 



the more accurate determination of the sample complete- 
ness were not behind this discrepancies, one possible expla- 
nation could come fr om the so-called evolution correction 
IIBlanton et al . 2003b), which we have not taken into account 
in the determination of our absolute magnitudes. Because we 
are dealing with relatively nearby objects, it seems unlikely 
that this correction could modify our results significantly. 
It is worth mentioning at this point that relatively small 
variations in the shape of the LF, which are probably not 
physically significant - given the uncertainties we are deal- 
ing with - , may translate into considerable changes in the 
values of the best-fit Schechter parameters. It is not conve- 
nient, therefore, to make comparisons between different LFs 
by just looking at these best-fit Schechter parameters. 

We have seen that the bright-end bump that we see in 
the blue bands (and marginally in the rest of the bands) is 
statistically very significant, according to our standard boot- 
strapping error analysis. Moreover, as we mentioned above, 
we also find c l ear - b ut noisy - evidence of its existence in 
iBlanton et aLI (|2003bh . We have also checked that this excess 
is not a consequence of any of the limits that we have im- 
posed to define our samples. This is, therefore, a remarkable 
result that may have strong implications for galaxy evolu- 
tion. 

From a preliminary analysis on the nature of this bright- 
end bump in the 01 u-band LF, we have seen that it is mostly 
populated by star- forming and active galaxies (~ 85%). The 




Figure 10. The 01 r-band SDSS DR6 Luminosity Function for 
blue and red galaxies separately. The SWML LF estimates are 
shown in diamonds. The dashed lines represents the best-fit 
Schechter function. Best-fit values of Schechter parameters a, M, 
and <E>* for both blue and red galaxies are also shown in the fig- 
ure. Shaded regions represent the la uncertainty calculated using 
a bootstrapping technique. 

spectra of these galaxies seem consistent with what we ex- 
pect from QSOs/Seyferts I (~ 60%) and LINERs/Starbursts 
(~ 25%). It seems, therefore, that an important fraction 
of the light that we receive from the brightest galaxies in 
the 01 ti-band would come from nuclear activity. Only about 
15% of galaxies in bright-end bump seem to be bulge-like 
galaxies. A more detailed study is needed, however, to fully 
understand the nature of this bright-end bump. 

The implications of this new results could be investi- 
gated using semi-analytic models of galaxy formation and 
evolution (SAMs). With these models, we can, in princi- 
ple, evaluate different processes and feed-back relations that 
could reproduce our results. 



5 SUMMARY 

In this work we make use of the SDSS Sixth Data Release to 
estimate the number counts, luminosity functions and lumi- 
nosity densities of galaxies in all SDSS photometric bands. 
The SDSS DR6 is, by far, the most complete survey of the 
nearby universe, containing redshifts for ~ 1, 000, 000 galax- 
ies down to magnitude r ~ 17.77 and covering ~ 7400 deg 2 
on the sky. The huge increase in the galaxy statistics with 
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respect to the previous work of iBlanton et all (|2003bh and 
the adequate treatment of both angular and brightness- 
dependent redshift completeness in our samples have allowed 
us to estimate the galaxy LFs of the nearby universe with 
unprecedented accuracy. In addition, we have calculated, for 
the first time, the SDSS galaxy number counts in all photo- 
metric bands using spectroscopic data. Luminosity densities 
in all SDSS bands have also been computed. 

The main results of this work can be summarized as 
follows: 



• The SDSS DR6 galaxy number counts in all SDSS 
photometric bands are consistent with an Euclidean, not- 
evolved Universe within a magnitude range that is limited by 
redshift incompleteness in the bright end and by the intrinsic 
apparent magnitude r-band limit of the survey, r = 17.77, 
in the faint end. 

• The SDSS DR6 LF of galaxies in the very blue 01 u- 
band d eviates considerably from that of IBlanton et al.l 
(2003b). This discrepancy can be explained by their lack of 
statistics, which has increased by a factor of ~ 9 thanks to 
the SDSS DR6. In the faint end, where we can reach about 1 
to 1.5 magnitudes deeper without losing accu racy, our SPSS 
DR6LF is considerably steeper than that of IBlanton et al.l 
(2003b). More interesting, in the bright end we find a re- 
markable excess, of ~ 1.7 dex at Mo.i u ~ —20.5 with respect 
to the best-fit Schechter LF. This bright-end bump is very 
strong in the 01 w-band and weakens in the 01 p-band, fading 
away towards the very red 01 z-band. 

• We conclude that the SDSS DR6 LFs of galaxies in 
the A q, 01 r, 0A i and °' 1 z bands are also compatible with 
IBlanton et ail l|2003bl ). considering the large increase in the 
statistics thanks to the SDSS DR6, of a factor ~ 4 — 5 in 
these bands. Some significant differences exist, however, es- 
pecially in the faint end, where we find a slightly steeper 
slope and we can reach about 1 magnitude deeper without 
losing precision. 

• A preliminary analysis of the origin of the bright-end 
bump seen in the °Vband SDSS DR6 LF reveals that it 
is comprised of QSO and Seyferts I galaxies (~ 60%), Star- 
burst and LINERs (~ 25%) and bulge-like galaxies (~ 15%). 
It seems, therefore, that a big fraction of this exceeding lu- 
minosity might comes from nuclear activity. 

• The Vband SDSS DR6 LF of blue galaxies is consis- 
tent with a Schechter LF with a remarkably steep faint-end 
slope, corresponding to a = —1.36. The 01 r-band SDSS 
DR6 LF of red galaxies has, however, a slightly decreasing 
faint-end slope, corresponding to a = —0.81. 

• The SDSS DR6 lumi nosity densities of gal axies are in 
very good agreement with IBlanton et all l|2003bf ) in all pho- 
tometric bands, since they are integrated quantities. 



The state-of-the-art results presented in this paper may 
be used to constrain a variety of aspects of star formation 
histories or feed-back processes in galaxy formation models. 
However, much effort is still needed in the survey field to 
fully understand the mechanisms that drive the evolution 
of galaxies in the Universe. This is especially necessary at 
high-z, where statistics are still very poor. 
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